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SOLAR  ATMOSPHERIC  DYNAMICS 

I .  INTRODUCTION 

The  goal  of  this  research  is  to  study  solar  atmospheric  waves  and  their 
role  in  heating  the  chromosphere  and  corona.  To  this  end  we  are  calculating 
the  propagation  of  radiating  acoustic  shock  waves  and  their  effect  on  solar 
spectral  lines.  The  ability  of  acoustic  waves  to  heat  the  corona  has  been 
cast  into  doubt  because  they  are  not  observed  to  carry  sufficient  energy  into 
the  transition  region  and  because  the  corona  is  observed  to  consist  of  magnetic 
flux  tubes.  To  clarify  this  issue  we  have  made  detailed  theoretical  studies 
of  acoustic  wave  flux  and  heating  of  the  chromosphere  and  are  developing  a 
"modal"  code  to  calculate  magnetic  type  wave  propagation.  As  a  preliminary 
step  we  have  studied  the  wave  vector  surface  of  Magneto-Acoustic-Gravity  Waves. 
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ll.  RADIATIVE  SHOCK  DYNAMICS 


The  development  of  computer  codes  to  study  radiating  shocks  was  begun. 

These  programs  will  be  applied  to  study  the  heating  of  solar  and  stellar  chromospheres 
by  acoustic  waves  and  the  detailed  analysis  of  spectral  line  profiles  in  the 
presence  of  acoustic  waves.  The  first  stage  in  this  project  was  the  development 
of  a  grey-LTE-continuum  radiation  hydradynamic  code,  which  is  described  below. 

The  dynamics  is  determined  by  the  equations  of  mass  and  momentum  conservation. 
Conservation  of  Mass: 


3R  1  . 

8m  p 

where  m  is  the  mass  column  density  of  position  R. 
Conservation  of  Momentum: 


80 

8t 


3m 


(P  +  Q)  -  g, 


where  Q  is  an  artifical  viscosity  to  stabilize  the  numerical  calculations. 

Q  =("  '  ax' 


to 


,  8U.  2  2  8U  / 

a  if  17  "v0 


otherwise 

where  l  is  a  constant  with  the  dimensions  of  a  length. 
Velocity  U  and  position  R  are  related  by 
3  R 


3 1 


=  U. 


These  equations  are  solved  by  dividing  the  atmosphere  into  layers  of  given 

mass  column  density  Am.  =  m.  .  -m.;  and  transforming  the  differential 

3  ^  1  /  2  3  1  3 

equations  into  finite  difference  equations  on  this  Lagrangian  grid  (Christy, 

196^;  Richtmyer  and  Morton,  1967). 

n+1/2 


The  velocity  U, 


and  position  are  defined  at  the  zone  boundaries, 


with  the  velocity  at  the  half  integral  times  tn+^^  =  tn  +1/2  A  t  and  the 


position  at  the  integral  times  t.  All  the  other  variables  are  defined  at 


The  finite  difference  equations 


the  zone  centers  j+1/2  and  integral  times  tn. 
are  solved  explicitly  for  Ujn+1/2,  then  Rjn+1 ,  then  p  j*{/2,  and  finally 
Q  by  making  a  sweep  through  the  grid  and  solving  the  finite  difference 

equations  in  the  following  order: 

First  solve 


n+1/2  _  u  n 

_J _ i_ 


1/2 


At 


n 


=  -g- 


pn  pn  „n-l/2  Qn-1/2 

“ j+1/2  -  j-1/2  *  ^ j+1/2  wj-l/2 

Am  . 


for  the  velocity  Ujn+1/2>  j  =  2 . .  N-l,  where  Am.  =  1/2  (  Amj+1/2 


Amj-l/2) 


,  n  .  /  . . n+1/2  ,.n-l/2x  ..  ..n+1/2  .n+l  .n 

and  Am.  ,  =  m .  ,  -ra .  and  At  =  1/2  (At  +  At  )  and  At  =  t  -  t 

j+l/2  J+1  J 

Second  solve 


Rn+1  -  Rn 


At 


- J-  =  U . 

n+1/2  j 


n+1/2 


n+1 


for  the  position  of  the  zone  boundaries  ,  j  =  1,  . .  N.  Third  solve 


R.?+1-  Rn+1 

_J±i _ .  i_ 

Am .  ,  P  n+1 


j+1/2 


j+1/2 


n+1  . 


for  the  density  p  j  =  1,  . .  N-l.  Finally  calculate  the  new  artificial 


viscosity 

Q  0+1/2 
y j+1/2 


=  ■  -2  (p  jii/2  ♦ p  j+1/2^  -  ”jn+1/2>2’ 


.  „  n+1  ^  n 

when  P.+1/2  >  PJ+1/2 

=  0  otherwise. 

We  take  2.  The  velocity  must  be  specified  at  the  two  boundaries.  We 

use  a  piston  at  the  lower  boundary 

U^n+l/2  =  Asin  (  wtn+1/2), 
and  a  transmitting  upper  boundary 

nn+l/2  n-1/2  ,„n-l/2  [,n-l/2v  ~n  A  tn  /  ( Rn  -  R  n) 

=  UN  (  N-l  UN  '  5N-l/2  1  N-i 


where  S  =  V  YP/P  is  the  sound  speed. 

The  thermodynamics  is  determined  by  the  equations  of  energy  conservation, 
radiative  transfer,  and  LTE  ionization. 

Energy  conservation: 


TT  +  (P+Q)  It  =  4lT^  (J_B)  " 
at  at  p  v  am 

where  E  =  3/2  nkT/  p  +  x  X  ^/m  is  the  thermal  plus  ionization  energy  per 

unit  mass,  and  n  is  the  total  particle  density,  x  is  the  fractional  ionization, 

X  is  the  ionization  potential  and  m  is  the  mean  mass  per  particle.  The 

conductive  flux  is  F  =  7.6  x  10  ^  T"3^2  ill 

Transfer  equation: 


2  a 


3  m 


"  1  ~|  3J 

k'+a  J  3  m 


(J-B) 


where  <  is  the  absorption  opacity  and  a  is  the  scattering  opacity,  and 
p2  -  1/3. 


Saha  equation: 

.«(!).*  m^'2  p-1.-,41 

where  K  =  0.333368  for  hydrogen  and  fA  is  the  fraction  of  electrons  from  heavier 
elements  (assumed  fixed). 

These  equations  are  solved  in  finite  difference  form  on  the  same  grid 
as  the  dynamics,  implicitly  and  simultaneously  be  Newton-Raphson  itleration 
(complete  linearization).  The  finite  difference  equations  are: 

Energy: 


E 


j 


n+1 

+1/2 


H  . 
J 


—  F  ^ 
J+l/2 


n+1 

+1/2 


A  t 


+ 


n+1 

j+l/2 


+ 


Q 


n+l/2l 

j+l/2  J 


=  o, 


n+1  n 

'j+l/2  p j+l/2 


n+1/2 
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where  the  heat  input  is  H  .  ,  =  4ir , 

J+1/2 


j+1/2  (Jj+l/2"Bj+l/2)  '  Fj^--  Fj-, 

j+1/2 


and  the  conductive  flux  is 


p j+1/2  +  Pj-l/2 


Fj  =  ‘  C,1/2  ( X j  +  1/2  +  X j-1/2  ^  Am.  iy,  +  Am, 


j+1/2 


T,i+l/2  *  T j-1/2 .  5/2 

O  ' 


( T j+1/2  '  T j-1/2 


j-1/2 

). 


The  boundary  conditions  on  the  flux  are  that  C  =  0  below  the  temperature  minimum 

and  FN  :  0  or  tn_3/2  =  TN-l/2  at  the  uPPer  Pounciary- 
The  number  density  of  particles  is 

n  =  n„  { l+x+B+A( 1+f )} 

H 

=  p  { l+x+B+A( 1+f )}  /m^  (l+4B+w) 

=  p  { l+x+B+A(l+f)}  /m 

where  B  is  the  ratio  of  helium  to  hydrogen  by  number,  A  is  the  ratio  of  heavier 
elements  to  hydrogen  by  number,  w  is  the  mean  atomic  weight  of  the  heavy  elements 
times  A,  and  f  is  the  fraction  of  electrons  from  heavy  elements.  The  total 
internal  energy  per  unit  mass  is 

E  =  Eth  +  Eion 

Eth  =  3/2  S'  T  (  l+x+B+A(l+f)} 

E.  =  =  T .  x,  T.  =  -£■  =  1.58  x  105  °K 
ion  m  ion  ’  ion  k 

Transfer: 


J.1+3/2  ~  J j+1/2  -  Jj+l/2  ~  ^ j-1/2 

ATj 


At 


j+1/2 


At  .  . 
J+l 


where  At. 

J+1/2 


-  X 


j+1/2  Amj+l/2 


^ J j+1/2  "  Bj+l/2}  =  ° 


At. 

J 


1/2  (  At  j+1/2  + 


ATj-l/2} 


JJ?  =  1/3 

j/[  =  p  »t)  f'rom  tabie 


B 
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The  upper  boundary  condition  on  the  radiation  field  is  no  incident  radiation, 
so 


JN-l/2  '  JN-3/2 


Tf 


N-l 


=  -1/2  (jn-1/2  +  JN-3/2)  +  BN-l/2 


At 


N-l/2 ‘ 


The  lower  boundary  condition  is 


J3/2  =  B3/2 

that  is  an  adiabatic  first  layer. 

Saha : 

Xn+*  -  2  <i>/  U  +Af+  (u<P+Af)2  +  i^}  =0 

J+l/2 

where  41  =  K  (k/m)  p  ^  T exp(  x^kT) 
and  K  =  0.333368,  X /k  =  1.58  x  105  °K. 

These  three  equations  are  solved  by  linearizing  them  in  the  corrections 
AT,  5J,  and  Ax,  and  solving  the  resulting  block  tri-diagonal  system  for 
these  corrections. 

After  the  new  temperature  distribution  is  obtained,  the  Saha  and  transfer 
equations  are  solved  for  the  ionization  x  and  mean  intensity  J  consistent 
with  that  temperature  distribution.  The  formal  solution  of  the  transfer  equation 
involves  sums  of  large  terms  (  At  )-i  and  small  terms  (  At  )  at  small  optical 
depth  and  in  the  usual  tri-diagonal  Gaussian  elimination  scheme  all  significance 
would  be  lost.  Therefore  to  be  able  to  carry  the  calculations  out  to  optical 
depths  as  small  as  the  smallest  machine  floating  point  number  the  elimination 
scheme  was  modified  as  follows:  The  transfer  equation  has  the  form 


Ai  Ji-1  +  Bi  Ji  +  Ci  Ji+1  =  Di. 

The  usual  procedure  is  to  calculate  elimination  matrices 

Ei.i  *  <Bi  *  *i  M  Ei.i  ■  -(Bi  *  *i  Ei-irl  <°i  *  «i  Ei-i>. 


where  these  are  then  used  to  back  substitute  for  the 


J. 

l 


E.  J.  ,  +  F. . 

l  l+l  l 


-9- 


At  small  optical  depth  the  radiation  field  is  nearly  constant  and  Ei  nearly 
the  identity  matrix  1.  It  is  therefore  possible  to  gain  accuracy  by  working 
with  a  matrix  c  i  equal  to  the  difference  between  and  1. 

That  is,  let 

ei  ■ 1  -  Ei 

Ji  ■  O  -  €  i>  Jt*i  *  Fi 

The  relations  to  calculate  and  are  found  by  substituting  this  equation 
into  the  difference  equation.  The  result  is 

ci  *  -  4i  'i-iV1  (Bi  *  *  ci  -  4i  S-i) 

Fi  =  (Bi  *  *i  -  H  ei.J  (,Di  -  Ai 
The  large  terms  in  the  sum  are  cancelled  analytically . 


III.  WAVES  IN  THE  CHROMOSPHERE  AND  TRANSITION  REGION 


Acoustic  pulses  are  seen  in  the  transition  region  ( AFGL-TR-77-Q108) , 
but  observations  of  line  widths  seem  to  rule  out  sufficient  acoustic  flux 
to  heat  the  upper  chromosphere  and  corona  (Athay  and  White,  1978;  Brunner 
and  McWhirter,  1979).  Dr.  Jorge  Vernazza  and  I  have  studied  the  propagation 
of  acoustic  waves  through  the  solar  atmosphere,  and  we  find  that  indeed,  acoustic 
waves  with  velocity  amplitudes  consistent  with  the  observed  wave  velocity 
of  1  Km/s  at  the  temperature  minimum  can  not  carry  sufficient  energy  to 
the  transition  region  and  the  corona  to  heat  them. 

We  have  used  a  one-dimensional  Lagrangian  fluid  dynamic  code,  including 

both  electron  conduction  and  optically  thin  radiative  losses  (a  la  Cox  and 

Tucker,  1969)  in  the  energy  equation.  The  waves  were  driven  by  a  sinusoidal 

piston  in  the  convection  zone,  1.5Mm  below  t  =  1,  and  propagated  through 

buuu 

an  empirical  solar  atmospheric  model  (Vernazza,  Averett,  Loeser,  1973)  - 

The  most  surprising  result  we  found  was  a  resonance  at  195s  period,  where 
the  velocity  amplitude  of  the  waves  at  the  temperature  minimum  for  a  constant 
driving  piston  amplitude  became  very  large  (Figure  1).  This  resonance  occurs 
near  the  acoustic  cutoff  frequency  appropriate  to  the  temperature  minimum, 
but  actually  seems  to  have  a  more  global  structure  and  may  be  the  chromospheric 
mode  found  by  Ando  and  Osaki  (1975).  This  "200s",  "chromospheric"  mode  extends 
higher  up  into  the  atmosphere  than  longer  period  modes  and  is  responsible 
for  the  change  over  in  the  dominant  oscillation  frequency  from  about  300s 
to  about  200s  as  one  goes  up  in  the  chromosphere  (Figure  2).  Both  the  driving 
frequency  mode  and  the  "200s"  mode  are  excited,  but  the  amplitude  of  the  long 
period  oscillation  decreases  rapidly  with  height  above  '£'  =  1,  while  the  amplitude 
of  the  "200s"  mode  is  nearly  constant  with  height  until  the  top  of  the  chromosphere. 
Oscillations  with  periods  greater  than  200s  excite  this  "chromospheric"  mode, 
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so  in  the  upper  chromosphere  only  this  "200s"  mode  is  observed.  Oscillations 
with  periods  shorter  than  200s  do  not  excite  the  "chromospheric"  mode  and 
have  large  amplitude  in  the  upper  chromosphere. 

The  implications  for  coronal  heating  of  our  results,  are  that  acoustic 

waves  of  all  periods,  50s  <  P  < 300s  with  velocity  amplitude  £  1  km/s  at 

ij  2 

the  temperature  minimum  carry  £  2x10  ergs/cm  s  into  the  transition  region. 

This  is  an  order  of  magnitude  less  than  the  required  energy  input  into  the 
corona  (Figure  3).  Even  increasing  the  driving  amplitude  does  not  significantly 
increase  the  acoustic  flux  in  the  transition  region.  Rather  larger  driving 
amplitudes  lead  to  greater  wave  dissipation  in  the  upper  chromosphere  (Figure 
4). 

Hence,  on  both  observational  and  theoretical  grounds  acoustic  waves  are 
ruled  out  as  the  energy  supply  for  the  corona.  The  heating  of  the  upper1 
chromosphere  is  still  an  open  question.  Magnetic  fields  clearly  play  a  role 
in  the  structure  of  the  upper  chromosphere.  But,  the  nature  of  that  role 
is  as  yet  unclear.  Magnetic  fields  may  alter  the  wave  generation,  tney  certair.lv 
change  the  types  of  waves  that  can  propagate  and  their  propagation  ar.d  ■■oi^r,:on 
properties,  finally  the  field  channels  heat  conducted  back  down  from  1  he  corona 
to  about  20.000K.  On  the  observational  side,  the  large  Ca  II  and  Mg  II  emission 
is  confined  to  the  regions  of  strong  magnetic  field  in  the  network.  However, 

Ca  II  H  and  K  bright  points  on  the  blue  side  of  the  line  are  fairly  uniformly 
distributed  over  the  solar  disk  throughout  the  interiors  of  cells,  and  are 
likely  due  to  upward  propagating  acoustic  waves.  As  one  step  in  exploring 
the  effects  of  magnetic  fields  we  have  calculated  the  shapes  of  wave-vector 
surfaces  for  magneto-acoustic-gravity  waves. 
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IV.  MAGNETO-ACOUSTIC-GRAVITY  WAVES 

The  solar  corona  and  chromosphere  require  some  non-thermal  energy  input 
to  heat  them.  This  energy  must  be  transported  by  motions:  either  wave  or 
quasi-static.  These  motions  are  produced  either  directly  by  the  turbulent 
convective  velocity  field  or  by  some  type  of  thermal  overstability.  Typical 
times  scales  are  comparable  to  or  longer  than  the  natural  buoyancy  time  scale, 
so  that  gravity  significantly  influences  the  motions.  Because  the  atmosphere 
is  compressible,  gas  pressure  influences  the  motions  as  well.  We  also  know 
that  the  structure  of  the  corona  is  controlled  by  the  magnetic  field,  which 
must  therefore  play  a  crucial  role  in  the  process  of  energy  transfer  and 
dissipation.  Dissipation  likely  takes  the  form  of  current  dissipation 
in  the  ionized  corona,  but  the  energy  that  drives  these  currents  must  be 
transferred  from  the  convective  zone  to  the  upper  atmosphere  by  wave  and 
quasi-static  motions.  (Except  for  currents  carried  up  into  the  corona  by 
emerging  flux,  which  is  a  transient  phenomena.)  Hence  the  energy  transport 
process  will  be  significantly  influenced  by  three  forces:  pressure,  gravity 

and  magnetic.  Nelson  Hartunian  and  I  explored  some  of  the  propagation  properties 
waves  under  the  influence  of  these  three  restoring  forces:  magneto-acoustic- 
gravity  (M.A.G.)  waves.  We  know  that  acoustic  waves  are  incapable  of  supplying 
sufficient  energy  to  heat  the  corona.  Can  the  broader  class  of  MAG  waves  supply  all 
the  non-radiative  energy  input  needed  to  heat  the  corona? 

Before  venturing  into  calculations  with  the  full  MAG  wave  equations, 
an  understanding  of  the  wave  properties  in  the  WKB  limit  is  a  helpful  guide. 

We  have  restricted  our  attention  to  the  case  of  constant  sound  speed  s  and 
Alfven  speed  a.  This  is  unrealistic.  While  an  atmosphere  can 
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have  constant  magnetic  field  §o  and  temperature  T,  an  atmosphere  in  hydrostatic 
equilibrium  must  have  density  p  decreasing  outward,  so  the  Alfven  speed 
a  =  B  /  (4up)  will  increase  outward.  By  neglecting  the  variation  of  Alfven 
speed,  we  restrict  ourselves  to  the  local  dispersion  relation  (McLellan  and 
Winterberg,  1968) .  This  is  valid  only  for  wavelengths  small  compared  to  the 
density  scale  height.  On  the  Sun  we  are  interested  in  waves  with  wavelengths 

comparable  to  the  scale  height,  so  we  must  integrate  the  full  wave  equation. 

However,  the  local  dispersion  relation  is  already  complicated,  and  a  picture 

of  its  properties  will  be  a  useful  guide  in  the  full  calculations.  The 

variation  in  Alfven  speed  with  height  was  considered  by  Yeh  (1974),  but  he 

only  examined  the  special  cases  of  horizontal  and  vertical  magnetic  field. 

The  local  dispersion  relation  has  been  studied  for  special  cases  by  several 

authors.  Bel  and  Mein  (1971)  and  Bel  and  Leroy  (1977)  have  discussed  the  case 

of  vertical  propagation  and  derived  the  variation  of  the  acoustic  cutoff 

frequency  with  magnetic  field  strength  and  direction.  The  special  cases  of 

propagation  either  in  or  perpendicular  to  the  g-B  plane  and  vertical  or 

-  ~o 

horizontal  B^  have  been  discussed  by  Chen  and  Lykoudis  (1972,  Michalitsanos 
(1973)  and  Nakagawa,  Priest  and  Wellck  (1973).  We  consider  the  general 
case  and  study  the  wave  vector  surfaces  of  the  local  dispersion  relation 
for  magneto-acoustic-gravity  (MAG)  waves  to  gain  insight  into  the 
behavior  of  these  waves.  The  wave  vector  surface  is  the 

three  dimensional  locus  of  the  real  part  of  the  wave  vector  for  all  waves 
at  a  given  frequency.  The  wavelength  is  given  by  |Re  k  |  =  2tt/X,  so,  for 
example,  small  wave  vectors  correspond  to  large  wavelengths,  and  since  the 
phase  speed  is  c  =  u/ |  Re  k  |  ,  waves  with  small  wave  vectors  propagate 
rapidly.  Each  type  of  restoring  force  tends  to  give  the  wave  vector  surface 
a  distinctive  shape  as  we  discuss  below.  We  illustrate  and  discuss  the 
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wave  vector  surfaces  of  MAG  waves  in  general  cases  in  the  next  section. 
Analytic  results  for  limiting  cases  of  (i)  high  frequency,  (ii)  vertical 
propagation,  (iii)  large  horizontal  wave  vector  kx  »  and  (iv)  horizontal 
magnetic  field  are  derived  in  an  appendix.  The  principle  results  are: 

(1)  At  high  frequencies,  where  buoyancy  is  negligible,  MAG  waves  behave 
like  magnetohydrodynamic  (MHD)  waves.  (2)  At  low  frequencies  in  a 
weak  magnetic  field  there  are  both  a  gravity  fast  mode  and  a  magnetic  slow 
mode,  but  these  exist  only  for  small  wave  vectors.  At  large  |kj  both 
their  wave  vector  surfaces  degenerate  into  planes  through  the  origin 
in  k-space  perpendicular  to  the  magnetic  field.  (3)  At.  low  frequencies 
in  a  strong  magnetic  field  there  is  a  magnetic  pressure  last  mode  with 
spheroidal  wave  vector  surface  about  the  origin  in  k-space.  The  gravity 
mode  is  suppressed  except  for  directions  nearly  perpendicular  to  the 
B  -  g  plane  or  for  nearlv  vertical  magnetic  field.  Awav  front  the  origin 

'O' 

both  modes’  wave  vector  surfaces  become  planes  through  the  origin  in  k-spa< 
perpendicu Lar  to  the  magnetic  field.  (4)  There  is  a  critical  frequency 
separating  high  from  low  frequency  behavior,  which  results  from  the  density 
stratification  of  the  atmosphere.  The  horizontal  component  of  the  magnetic 
field  increases  the  density  scale  height  and  thus  reduces  the  critical 


f  requency . 
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A<  WAVE  VECTOR  SURFACES 


The  local  dispersion  relation  for  magneto-acoustic-gravity  (MAG) 


(McLellan  and  Winterberg,  1968).  Here  B^  is  the  unperturbed  magnetic 
field,  is  the  ambient  density,  is  the  horizontal  component  of  the 

wave  vector  k,  s  =  (yp  /p  )  2  is  the  sound  speed,  a  =  B  /(4irp)  is  the 

-  0  0  .00 

Alfven  speed,  and  y  is  the  ratio  of  specific  heats.  We  solve  this  for 

k  (to,  g,  B^) .  The  first  factor  is  quadratic  in  k  and  represents  the  usual 

Alfven  waves,  which  we  discuss  no  further.  The  second  factor  is  quartic 
in  k  and  represents  the  MAG  waves.  In  the  limit  of  zero  magnetic  field  it 
reduces  to  the  dispersion  relation  for  acoustic-gravity  waves.  In  the  limit 
of  zero  gravity  it  reduces  to  the  dispersion  relation  for  magneto-hydrodynamic 
(MHD)  waves,  and  in  the  limit  of  zero  compressibility  (y,  s,  H  =  PQ/p0g  ~>c°) 
it  reduces  to  the  dispersion  relation  for  magneto-gravity  waves. 

The  wave  vector  surface  is  the  three-dimensional  surface  of  those 

values  of  the  real  part  of  the  wave  vector,  Re  k  (w,  g,  B^) ,  that  satisfy  the 

dispersion  relation  for  a  given  frequency  u.  Energy  propagates  perpendicular 

to  this  surface.  We  will  study  three-dimensional  pictures  of  the  wave 

vector  surface  and  two-dimensional  slices  through  the  surface  to  help 

clarify  its  shape.  One  can  recognize  the  type  of  wave  mode  from  the  shape 

of  its  wave  vector  surface.  Where  pressure  is  the  dominant  restoring  force 

(acoustic  or  magnetic  pressure  modes),  there  is  no  preferred  direction  and 

the  wave  vector  surface  is  a  spheroid  about  the  origin  in  k-space. 
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Where  buoyancy  is  the  dominant  restoring  force  the  vertical  (-g)  is  the 
preferred  direction.  Gravity  produces  a  characteristic  acceleration  for 
each  direction  of  fluid  displacement.  Each  direction  of  motion  therefore 
corresponds  to  a  particular  frequency.  Because  the  fluid  is  not  compressed 
the  motion  is  transverse,  k.  u  =  0.  Since  at  any  given  frequency  motion  is 
possible  only  at  a  single  angle  to  the  vertical,  the  wave  vector  must  make 
the  complementary  angle  to  the  vertical  and  the  wave  vector  surface  is  a  cone 
of  revolution  about  the  vertical  axis.  Where  magnetic  tension  is  the  dominant 
restoring  force  the  magnetic  field  ,  is  the  preferred  direction.  Magnetic 
field  lines  act  like  rubber  bands  under  tension.  Waves  are  transverse  and 
propagate  along  the  field  lines  with  a  characteristic  speed  equal  to  the 
square  root  of  the  tension  divided  by  the  density.  The  wave  vector  surfaces 
are  planes  perpendicular  to  B^  • 

The  character  of  the  MAG  waves  depends  on  their  frequency  and  the 

magnetic  field  strength.  Frequency  is  scaled  by  N  =  yg/2s,  the  acoustic 

cutoff  frequency  in  the  absence  of  a  magnetic  field.  For  high  frequencies, 

w/N  >>1,  gravity  is  unimportant  because  the  fluid  acceleration  is  much 

ac 

larger  than  g,  so  the  waves  behave  as  fast  and  slow  mode  MHD  waves  (see  appendix 
eqns  A4,  A5).  Acoustic  waves  cannot  propagate  at  frequencies  below  the  acoustic 
cutoff  frequency.  For  low  frequency  waves,  w/N  <1,  gravity  is  significant. 

3.C 

Magnetic  field  strength  is  parametrized  by  the  ratio  of  Alfven  to  sound  speeds, 
hi  ^ 

a/s  =  (B  /4ryp  )  -  (P  /P  ,  .  For  weak  fields  we  expect  to  find  an 

o  o  mag  gas; 

acoustic-gravity  mode  plus  a  slow  magnetic  mode.  For  a  strong  field,  we 
expect  to  find  MHD  waves  in  the  high  frequency  limit  and  a  new  situation 
where  all  three  restoring  forces  are  important  in  the  low  frequency  limit. 

We  now  consider  these  cases,  starting  with  the  simplest  situation  -  high 
frequencies,  where  gravity  is  unimportant.  Detailed  analytic  results  are 
derived,  where  possible,  in  the  appendix. 

1 .  High  Frequency 

High  frequency  means  frequency  above  the  acoustic  cutoff  frequency,  where 
gravity  is  negligible.  In  the  absence  of  a  magnetic  field  the  acoustic  cutoff  is 

N  =  yg/2s  =  s/2H,  (2) 

aC 

where  H  =  ?o/p^g  is  the  pressure  and  density  scale  height.  A  magnetic  field 
reduces  the  acoustic  cutoff  frequency  because  it  increases  the  scale  height. 

(A  cutoff  marks  a  transition  between  waves  that  propagate  vertically  and  those 
that  don't,  but  instead  vary  exponentially  with  height.) 
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Three  dimensional  pictures  of  the  wave  vector  surfaces  for  frequencies 

above  the  acoustic  cutoff  frequency  are  shown  in  figures  5  and  6  for  the 

cases  of  strong  and  weak  magnetic  fields  respectively.  Cuts  through  the  wave 

vector  surfaces  at  45°  to  the  B^  -g  plane  are  shown  in  figures  7  and  8.  In 

both  cases,  as  expected,  the  surfaces  look  like  the  usual  MHD  surfaces 

because  for  N  the  effect  of  gravity  is  negligible.  The  spheroidal 
ac  ' 

wave  vector  surface  near  the  origin  is  the  fast  mode  and  the  planar  wave 
vector  surfaces  on  either  side  of  it  are  the  slow  mode.  In  a  weak 
magnetic  field:  the  fast  mode  is  an  acoustic  mode  and  has  a  cutoff 
frequency,  below  which  it  ceases  to  propagate.  The  slow  mode  is  a  magnetic 
tension  mode  and  propagates  for  all  frequencies,  so  it  has  no  cutoff 
frequency.  In  a  strong  magnetic  field:  the  fast  mode  is  a  magnetic 
pressure  mode  and  propagates  for  all  frequencies,  so  it  has  no  cutoff 
frequency.  The  slow  mode  is  a  one-dimensional  acoustic  wave  propagating 
along  the  magnetic  flux  tubes  and  has  a  cutoff  frequency,  below  which  it 
ceases  to  propagate. 

The  acoustic  cutoff  frequency  can  be  found  by  considering  the  special 
case  of  vertical  propagation  (see  Appendix  and  Bel  and  Mein,  1971;  Bell  and 
Leroy,  1977).  In  the  weak  field  limit,  the  fast  mode  has  the  cutoff 
frequency  (eqn.  A 1 2 ) 

u>  =  N  s/(s^  +  a^sin^9) 2,  s>>a.  (3a) 

ac  ac 

Below  this  critical  frequency  the  fast  mode  becomes  an  internal  gravity  mode. 
In  the  strong  field  limit,  for  a  non-horizontal  magnetic  field,  the  slow 

mode  has  the  cutoff  frequency  (eqn.  A20) 

ijj  =  N  cos9,  a>>s,  9  <  tt  /  2 .  (3b) 

ac  ac 

Here  0  Is  the  p>lar  angle  of  the  magnetic  field  B^.  Below  this  critical 
frequency  the  slow  mode  does  not  propagate  except  in  a  few  special  directions. 
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The  minimum  value  of  the  acoustic  mode  cutoff  frequency  for  any  magnetic 

field  strength  occurs  for  a  horizontal  magnetic  field  and  is  (A22 ,  A43) 

oj  =  N  s/(s2  +  a2)^.  (3c) 

ac  ac 

These  result  on  the  cutoff  frequency  must  be  treated  with  caution,  because 
the  WKB  approximation  is  invalid  near  the  cutoff. 

The  acoustic  mode  has  im  k =  -i/2H  (fig  8b>  and  eqns.  A12 ,  A19,  A38) ,  as 

we  would  expect  from  our  experience  with  acoustic-gravity  waves.  The 

magnetic  mode  has  Im  -  0  (fig  8b,  and  eqns.  A14,  A16) .  These  are  general 

properties,  that  are  approximately  true  for  all  magnetic  field  strengths. 

They  are  a  consequence  of  energy  conservation.  The  wave  energy  flux, 

_  I  pu2s  (acoustic) 

F  =  pu  V  =  )  „  (4) 

group  ^pu  a  (magnetic'', 

must  be  constant  with  height.  For  acoustic  waves  in  an  isothermal 
atmosphere,  where  s  =  constant. 


u 


cc 


-z/2H. 

e 


(5) 


The  WKB  expression  for  the  velocity  amplitude  is 


-h  ik  z 


e  z 


and 


|  k  |  -  u/s  =  constant. 


Im  k  =  -  1/2H  (acoustic). 

z 


For  magnetic  waves. 


-h  ~  z/2H 
a  «  p  «  e  , 


u  a  p 


-1/A 


(6) 


(7) 

(8) 


(9) 

(10) 
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but 

g. 

j  k  j  -  w/a  p  2. 


(ID 


Hence  the  height  variation  of  the  velocity  amplitude  of  magnetic  waves  is 
due  to  the  variation  oi  1  k  '  “  and  not  the  exponential  factor,  so 

I  m  k  =  0  (magnetic).  (12) 


Different  modes  will  couple  with,  each  other  when  Mu  ir  frequencies  and 
wavelengths  are  the  same.  For  large  ;  k  1  there  is  at  most  one  upward  and 
townward  ;  ropagating  wove,  so  no  coupling  of’  soar,  bu‘  for  sa.al  !  1. 

coup ;  i n..;  is  sometimes  possible.  lit.-  magnetic  and  acoustic  modes  v :  1  .  have 
intersecting  wave  vector  surfaces  when  a>  ■  consider  tho  lit  it  in.'  c  .se 
of  ”eiticaliy  propagating  waves.  ;u  coupiin  :  condition  can  be  1 ound  by 
equating  tho  expressions  lor  the  col  parts  •  .  the  fast  and  slow  wave  vec  or 
(eqns.  Alb  and  A19) ,  and  solving  for  u>.  Resonant  coupling  occurs,  when 


a>>s,  for 


w/N  =  cosft  (1  +  s'  cos  G/a~ ) 
ac 


Re  k 

z 


cost  s 
2H  a  ’ 


(13) 


that  is  tor  w  slightly  above  the  cutoff  frequency  and  wavelength  larger 
than  the  scale  height.  Bel  and  Leroy  (1977 )  have  obtained  corresponding 


results  for  a  vertical  field  of  arbitrary  strength. 


2.  Low  Frequency 


We  now  turn  our  attention  to  the  more  interesting  but  complicated  case 
of  low  frequencies  where  all  three  forces  are  significant. 

(a)  Weak  Magnetic  Field 

In  a  weak  magnetic  field  one  mode  must  be  a  gravity  wave  and  the  other 
a  slow  mode  magnetic  tension  wave,  in  order  to  go  continuously  to  the  limit  of 
acoustic-gravity  waves  in  the  absence  of  a  magnetic  field  and  MHD  waves 
in  the  absence  of  gravity.  This  is  indeed  the  case,  as  shown  in  the  three- 
dimensional  wave  vector  surface  pictures  (fig.  9)  and  slices  through  the 
wave  vector  surface  in  figure  10.  There  are  some  modifications,  however.  First, 
the  gravity  mode  hyperbola  is  cut  off  near  the  Alfven  planes.  This 
is  similar  to  the  wave  vector  surfaces  of  magneto-gravity  waves  (Lighthill, 
1967;  Schwartz  and  Stein,  1975).  Second,  the  slow  mode  wave  vector  surface 
has  a  bulge  away  from  the  origin  in  the  vertical  direction,  rather  than  a 
dimple  toward  the  origin  in  the  magnetic  field  direction  as  found  in  MHD  waves. 
Third,  at  large  horizontal  wave  vector,  kj_  >>u)/min  (s,  a),  both  modes  wave 
vector  surfaces  become  planes  perpendicular  to  the  magnetic  field  direction 
whose  extensions  pass  through  the  origin.  This  can  be  seen  most  easily  from 
the  analytic  results  derived  in  the  appendix  (eqn.  A40)  and  from  slices  through 
the  wave  vector  surface  when  the  magnetic  field  is  nearly  vertical  (fig.  11). 
When  a  wave  vector  surface  becomes  a  plane  through  the  origin  perpendicular 
to  the  magnetic  field,  its  group  velocity  goes  to  infinity,  since  k  is 
independent  of  w,  and  its  potential  energy  becomes  much  larger  than  its 
kinetic  energy.  This  occurs  below  the  critical  frequency 

co  =  NcosG  (14) 

c 
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k  . 

(A36,  where  N  =  (y-1)  g/s  is  the  Brunt  Vaisala  frequency)  for  both  the  fast 
gravity  and  slow  magnetic  modes  at  sufficiently  large 

|  k  |  .  The  propagation  direction  perpendicular  to  the  Bq  -  g  plane  is  a 
special  case.  In  this  direction  the  plane  through  the  origin  in  k  -  space 
perpendicular  to  the  magnetic  field  has  zero  vertical  component.  Hence  in 
this  direction  at  large  kj_  the  wave  vector  surfaces  arch  over  until 
Re  k^  =  0  (fig  12).  In  this  special  case  there  is  a  sharp  cut-off  to 
propagation.  However,  in  general  whenever  the  wave  vector  surface  approaches 
a  plane  through  the  origin  perpendicular  to  the  magnetic  field  (Re  6  -*■  o  in 
A31)  the  waves  become  non-propagating. 


(b) Strong  Magnetic  field 

As  the  magnetic  field  strength  increases,  the  distance  (if  the  All'ven 

planes  from  the  origin,  u/a ,  decreases.  Any  gravity  mode  wave  vector 

surface  is  restricted  to  lie  between  these  planes,  on  which  magnetic  tension 

provides  all  the  needed  acceleration.  Hence  the  gravity  node  wave  vector 

surface  gets  squashed  (figs  13-15).  Instead  of  a  hyperboloid,  the  wave 

vector  surface  degenerates  into  two  intersecting  pianos  through  the  >riv.in  - 

one  perpendicular  to  the  magnetic  field  and  the  other  with  k  =  0.  file 

z 

critical  frequency  for  both  zero  and  large  kx  is  the  same  in  the  st  ron.. 

;  iel.i  limit  (eqn->.  A20  and  Aii) 

.  -  N  cos-.  .  no 

Hence  a  gravity  wave  slow  mode  cannot  propagate  in  a  strong  magnetic  f :  e*  id , 
except  for  a  small  range  of  intermediate  kj_  and  directions  marly 
perpendicular  the  B  -  g  plane  (fig  1 6 )  or  nearly  vertical  magnetic  field 
(fig  17).  These  last  two  figures  also  clearly  ilium  rate  the  degeneration 
cf  the  wave  vector  surfaces  at  large  into  planes  through  the  origin 
perpendicular  to  the  magnetic  field.  Again  Ke  k^  -  0  on  this  plant1  in  tilt 
special  direction  perpendicular  to  the  B  -  g  plane. 

The  wave  vector  surface  of  the  second  mode  -  the  magnetic  pressure 
fast  mode  -  remains  a  spheroid  about  the  origin  at  low  frequency  (i  ies  la- If 
as  at  high  frequency  (fig  5).  However,  at  low  frequency  this  spheroid 
attaches  onto  a  plane  through  the  origin  perpendicular  to  the  magnetic 
field  at  large  ,  rather  than  terminating  as  at  high  frequencies.  These 

waves  with  large  k_x  have  virtually  all  their  energy  in  the  form  of  magnetic 
potential  energy. 
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(c)  Conclusion 

The  significance  of  these  results  for  solar  heating  are:  First, 
all  low  frequency  waves  (u>  <  NcosQ)  in  a  strong  magnetic  field  propagate 
energy  in  the  direction  of  the  magnetic  field.  Thus  refraction  will  not 
reduce  the  energy  flux  of  low  frequency  waves.  Second,  the  group  velocity 
of  such  waves  may  be  greater  than  the  Alfven  speed.  And  third,  the 
magnetic  energy  density  will  be  greater  than  the  kinetic  energy  density. 
Because  of  these  last  two  results,  the  total  wave  flux  will  not  be 
observable  in  line  widths  or  Doppler  shifts.  Hence  current  observations 
do  not  place  any  restrictions  on  the  flux  of  low  frequency  M.A.G.  waves. 
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Appendix:  Limiting  Cases 


The  local  dispersion  relation  for  magneto-acoustic-gravitv  waves 
(neglecting  gradients  in  the  sound  and  Alfven  speeds)  is(McLellan  and 
Winterberg,  1968) 


4  2 

u>  -  u 


2  2  2 

(s  +  a  )  k  -  iyg  .  k 


2,2  (k  '  V  2  2 

S  *  "~47p -  +  (y_1)  ?  k- 


iyk 


(k  .  B  )  (g  .  B  ) 
-o  °  -o 


4  TT  P 


0, 


( A 1  ) 


where  B  is  the  ambient  magnetic  field,  p  is  the  ambient  density,  kx  is  the 
horizontal  component  of  the  wave  vector  k*?  and  y  is  the  ratio  of  specific 
heats.  This  dispersion  relation  is  a  fourth  order  equation  for  the  complex 
vertical  component  k  of  the  wave  vector  in  terms  of  the  frequency  and 
horizontal  component  kx ,  which  are  real: 

D  =  ak  ^  +  bk  ^  +  ck  ^  +  dk  +e=0  (A2) 

z  z  z  z 
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(A3) 


9  is  the  polar  angle  between  the  magnetic  field  and  the  vertical  and  $  is  the 
azimuthal  angle  between  the  horizontal  projections  of  tile  magnetic  field  and 
the  wave  vector-vertical  plane. 


(i)  High  Frequency 


For  high  frequencies,  u)>>g/s,  gravity  is  negligible  and  the  M.A.G.  waves 
become  the  well  known  MHD  waves.  The  dispersion  relation  reduces  to 

4  22  22  22  ~ 

oi  -  a)  (s  +  a  )  k  +  s  k  -■  ■  =  0  (A4) 

4'irp 

o 

or 

4  2  2  2  2  2  2 

c  -  c  (s  +  a  )  +  s  a  cos  0=0 


where  c  =o)/k  is  the  wave's  phase  velocity  and  0  is  the  angle  between  k 
and  B  .  The  wave  modes  ir.  this  case  are  the  usual  fast  and  slow  modes. 
Their  phase  velocity  is 


,  2  2  2 
(s  +  a  )  - 


,2  2  2. 

-4  s  a  cos  9 


kB 


■5 


(A5) 


Ln  the  limits  a»s  or  s>>a,  the  phase  velocities  are 


fast 


2  2.2 

c  +  c  s  m  o 


kB 


\s  Low 


2  2  , 

C<  C°S  9kB 
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where  c  =  max  (s,  a)  and  c,,  =  min  (s  ,  a)  (see  for  instance  Bazer  and 
Fleischman,  1959.)  In  the  weak  field  limit  the  fast  mode  is  an  acoustic 
mode  and  the  slow  mode  a  magnetic  tension  mode,  while  in  the  strong  field 
limit  the  fast  mode  is  a  magnetic  pressure  mode  and  the  slow  mode  a  one 
dimensional  acoustic  mode  propagating  along  the  magnetic  flux  tubes. 


(ii)  Vertical  Propagation 

For  vertical  propagation,  kj_  =  0,  the  dispersion  relation  reduces  to 
(Bell  and  Mein,  1971) 

Or/  -  Or)'  I?  -  s2  (1  +  4ff>J  <A7> 

2  2  i  2 

+  a  s  (1  +  T-jj-)  cos  0  =  0 


where  H  =  P/pg  =  s  /yg  is  the  pressure  scale  height  in  an  isothermal  atmosphere 


2  2 

This  is  similar  to  the  MHD  dispersion  relation,  with  s  s  (1  +  i/kH) .  The 


phase  velocities  are  now  a  function  of  k 


.2  a2  +  s2(l  +  i/kH) 
C  +  =  - 2 - 


t*  [<a2  +  s2)2  -  (fj)2  +  2i  ^ 


2  2 
(sZ  +  a) 


4  (1  +  ■—)  s2a2cos29 


(A8) 


In  the  weak  and  strong  field  limits  this  becomes 

2 


fast 


(*)♦ 


(1  +  ^  )  +  a  sin  9 


s>>a 


(A9) 


2  ^  2  ,  i  ,  .2. 

a  f  s  (1  +  — )  sin  0 


a>>s 


“slow 


=  fe)-2- 


a2s2(l  +  i/kH)cos29 


a2  +  s2  (1  +  i/kH) 


(A10) 


2  2 
a  cos  9 


s>>a 


2  i  2 

s  (1  +  — )  cos  0  a>>s 


These  can  be  solved  for  the  wave  vector  k.  Consider  four  cases  separately, 
(a)  Fast  mode  with  s>>a: 


2  9  2  2  9  2 

k  (s  +  a  sin  9)  +  i(k/H)  s  -  w  =  0  , 


(All) 


2H  2  2.2 

s  +  a  s  m  0 


2  2 
\Js  +a  s 


.  2 
sm  0 


1  - 


N  2 
ac 


2  2  2  .2 
oj  s  +a  sin  f 


(A12) 


where  N 


yg/2s  =  s/2H. 
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Tiure 


is  a  cutoff  frequency, 


crit 


N  s/ 
ae 


beli'w  which 


2  2 
s“  +  a^sin 


the 


wave 


is  evanescent. 


(AH) 


This  mode  is  a  cross  between  the  acoustic  branch  of  acoustic-gravity  waves 
and  the  fast  mode  Mill)  wave.  The  stratification  produce-;  the  cutoff 
frequency  and  a  vertical  amplitude  growth  (imaginary  part  of  the  wave 
vector i .  The  magnetic  field  reduces  the  cutoff  frequency  below  N 
and  reduces  the  vertical  growth  rate  below  !/2H,  their  values  for  "  a.-oustic- 
gravity  waves. 

<b>  Slow  mode  with  s*->a 


coso 

This  is  the  usual  magnetic  tension  slow  mode.  Gravity  does  not  alter  th-- 
n;oe»- . 


fast  trade  with  a>  s 

2  2 

in  the  case  where  kll^s^/a* 


2  ?  2  2  7  2  a 

k  (a"  +  s  sin^O)  +  i(k/H)  s  sin  b  -  or-  =  0 


s  o 


K+= 


2  2 

i  s  sin"1-- 
2H  2 


J  2  2  .  2 

va  +s  sin  ■ 


N  2 


w 


2 

s  sin  0 

2 

a 


(Alb) 


This  is  a  magnetic  pressure  mode  slightly  modified  hv  gravity.  There  is 
a  small  vertical  growth  in  amplitude  and  a  slight  dependence  of  the  phase 
velocity  on  the  direction  of  the  magnetic  field.  There  is  no  culjofjJ 
frequency  contrary  to  appearances,  because  as  cu  -*  (  ~j,  kH  v  s‘‘/a“  and 

our  initial  assumption  is  no  longer  valid.  The  case  o)  small  kH  must  be 
invest igated  separate  1 v . 


ii 
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2  2 

The  special  case  kH  <  s  /a  corresponds  to  small  frequencies,  m/N  <  s/a. 
In  this  case  the  vertical  dispersion  relation  reduces  to  ac 


Cl)3-‘  (f)  a 


+  i 


2  2  2 
s  a  cos  6 

mH 


=  0, 


which  can  be  written  as  a  cubic  equation  for  the  wave  vector 


k3  -  k 


u>4H 


2  2 
a  cos  9 


2  2  2« 
a  s  cos  0 


=  0  . 


(A17) 


In  the  limit  of  small  m,  the  first  and  second  terms  dominate,  and 


k  =  + 


acos0 


If  we  include  the  first  order  correction,  we  find  that  the  fast  mode  at  very 
low  frequency  in  the  strong  field  case  is 


k 


+ 


+ 


U) 

acos6 


+ 


(A18) 


The  third  root  of  the  cubic  is  pure  imaginary.  This  is  a  magnetic  tension 
mode  and  has  no  cutoff  frequency. 

(d)  Slow  mode  with  a>>s: 


i  2  2 

k  s  cos 


4-  i 


(k/H)  s^cos^0 


2 

0) 


=  0  , 


so 


k 


0) 

SCOS0 


1 


N  2  2. 

ac  cos  9 

2 

U) 


(A19) 


This  is  an  acoustic  mode.  Its  phase  velocity  in  the  direction  of  B  is  the 
sound  speed.  It  grows  in  the  vertical  direction  with  the  usual  sca?e  height  of 
2H.  It  has  a  cutoff  frequency 


U) 


crit 


N  cosO. 
ac 


(A20) 


The  case  of  a  horizontal  field  must  be  treated  separately.  For  0  =  90°  one  mode 
has  zero  phase  velocity  and  the  dispersion  relation  for  the  other  mode 
becomes 


,  2  ,  2.  2  ik  2 

(a  +  s)k  -  —  s  =  0 

H 


k  =  - 


i 

2H 


2  2 
a  +  s 


\n  2 

Na  +  s 


1  - 


2,2 

a  +  s 


(A21) 


In  this  case  the  cutoff  frequency  is 


a)  ,  =  N 

crit  ac 


2  2 
a  +  s 


(A22) 


and  the  vertical  growth  rate  is  reduced. 


(Hi)  Large  kj 


In  the  limit  of  large  kx,  that  is  kx>>  uj/min(a,s)  t£e  dominant 
terms  in  the  dispersion  relation  are  the  terms  of  order  k  : 


,4  22  2  3  .  2  2 

k^  as  cos  e  +  2kz  kxa  s  cosesm0cos() 


2  2  2  2  2  2  2 
+  k  kj_  a  s  (cos  9  +  sin  9cos  <j>) 


3  2  2 

+  2k^  kj_  a  s  cos9sin9cos<j> 


2  2  2  2  2 
=  (k  +  kj_  )  (k  cos9  +  kj_sin9cos()))  a  s 


2 ,2 (-  ■  y 

=  S  k  — : - 

4np 


0. 


(A23) 
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Henee  either 


or 


k  .  B  ->  0, 

-o 


2  2  2 
k  =0,  i.e.  k  =  -  k»  . 

-  z 


(A24) 


(A25) 


The  next  lower  terms  are  the  terms  of  order  k  ,  and  the  dispersion 
relation  is 


■) 

\C 


(k  .  b  r 

~o 


4  up 


o 


(k  .  B  )  (g  .  B  ) 

_j; _ ;q_ _ , _ o 

4  n  i !  2 

o  s 


„  k  .  B  r 

=  k  ——A1  k.  B  -  irA.li  /s 
4tio  ~o  o 


2  2 

=  (k  +  k_^  )  (k  cos6  +  kAsin'>cos-! ) 


k^_  cos0  +  k  sinOcoS't  +  icosB/H  a“ 


=  0  . 

The  four  roots  are 


(A If-  ) 


(1)  k  .  B  =  0 

~  ~o 


or 


k  =  -  k,  tan0cos<t> 


(A27) 


<2>  £  •  ?0  =  iy%  '  ?o/s 


or 


or 


k  =  -  kj_tan0cos<j>  -  i/H 


(3  and  4)  k^  =  0 


k  = 
z 


+  ik^ 


(A28) 


(A29) 
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The  first  two  modes  have  wave  vector  surfaces  that  are  planes  perpandicular 
to  g  ,  while  the  last  two  are  non~propagating2modes .  For  the  first  two 
modes  we  can  include  the  next  lower  order,  k  ,  terms  in  the  dispersion 
relation  as  a  perturbation.  The  dispersion  relation  now  becomes 


k  .  B  F  2 

k^  1 - Z°  k.  B  -  iYg  •  B  Is 

4irp  ~  -o  ~  ~o 


+  N2  kA2  -  k2  to2  =  0, 


wnere  N  =  (y  -  1 g/s  is  the  Brunt-Vaisala  the  frequency  of 

oscillation  of  a  fluid  parcel  in  an  isothermal  atmosphere. 

In  component  form 


(k  2  +  k,2)  (k  cos0  +  k.  sin9cos(J>) 

7.  ~  7.  •*- 


k  cos6  +  k.sin6cos<f>  +  icos6/H 

L2 


■] 


2  2  2 
,,  2  ,  2n  s  +  a  2  ,  N 

-  <kz  +  ki. )  — "  +  ~2 

s  a  a 


k,  =0. 


Let  us  define  6  by 


k  =  -  ki  tan0cos<j>  +  6. 
z 


then 


k  cosO  +  k.  sin9cos<j>  =  6cos8 

7. 


and 


k  2  +  ki2  =  k,2  ( tan20cos2<j>  +  1)  -  26k ,  tan0cos<j> . 

z  * 


The  dispersion  relation  becomes  (note  6~  H  ) 


6  cos20  (6  +  i/H)  -  w2  (a2  +  s2)  /a2s2 


+  (N2/a2)  /  (tan26cos2<{>  +  1)  =  0. 


(A30) 


(A31) 


(A32) 
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Hence 


6 


2 

u 


COS  0 


2  2 
s  +  a 

2  2 
a  s 


a  1  -  sin  0s in  $ 


(A3  3) 


There  exists  a  critical  frequency 
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2  2 
a  s 

2  2 
a  +s 


2 

COS  0 
4H2 


2  2 
1  -f  tan  9 cos  ' 


> 


2  2 
a  s 

2  2 
a  +s 


2 


a 


2 

cos 


0  • 


(A34) 


The  wave  vector  has  a  real  part,  <5cos9,  parallel  to  the  magnetic  field  for 

oj>oj  .  However,  for  td  <  id  ,  Re6  =  0  so  the  real  part  of  the  wave  vector  is 
c  c 

perpendicular  to  the  magnetic  field  through  the  origin. 


2  2  2 

w  =  N  cos  0,  (a>>s,  independent  of  .  (A351 

c  ac 

Note  this  is  the  same  as  the  cutoff  frequency  (A20)  of  the  slow  acoustic 
mode  in  the  opposite  limit,  k^K). 


For  a  weak  field. 


2 


Id 

C 


=  N2/  (1  +  tan*!)  cos2  4>) 
2  2 

>  N  cos  0. 


(s>>a) 


(A36) 


In  this  case  6  has  a  real  part  for  all  if  u)>N,  but  only  for  <j>  inside 
a  fan 

2  (  N2  A  ,  2 

cos  >  I  — 2  -  11 /tan  0 


if  N>id>Ncos0. 


(A37) 
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Above  the  critical  frequency 

(for  high  frequencies  or  nearly  horizontal  ,cos6-o) 


<5  =  -  “77  + 


2H  ± 


U) 

2  x  2 

s  +  a 

COS0  ' 

V  2  2 
s  a 

CO 

acos6 

s>>a 

0) 

SCOS0 

a»s 

MHD  slow  modes 

(A38) 


exponential  growth  factor. 


Below  the  critical  frequency  (uj<<coc),  with  a>>s 


+i/«  %  IT  - 30 

n  y  1  -  sin  Csin  ^ 


(A39) 


Hence  the  wave-vector  surfaces  are  planes  through  the  origin  and 
perpendicular  to  B 


Low  the  critical  frequency  (u>-;<  w  )  with  s>>a 


,  .  .  n  _ i _ _  i_ 

-  1  a  .  2„  2  h  2H  ’ 

(1  -  sin  Ssm  <J>) 


(A40) 


1  s  Y-l 
Ha  2 


Again  the  wave-vector  surfaces  are  planes  through  the  origin  perpendicular 


to  B  . 

~  o 
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(iv)  Horizontal  Magnetic  Field 

When  the  magnetic  field  is  horizontal,  the  dispersion  relation  reduces 
to  a  quadratic  in  (since  cos6  =  0) 


4  2  T,  2  ,  2,  J- 

o-o)  |_(  s  +  a  )  k 


2  2  2  2  2 
+  (N  +  aV)  sk/ 


-  iyg.k 


2.2  2.  2  2. 

o  (s  +  a  )  +  s  a  k 


2 

COS  (J> 


-  iyg  o)  k 


4  2.2  2 

4-  o)  -  o)  (s  +  a 


)„/  -  ,V  kA‘  +  »2Aa2] 


For  small  k^and  for  all  kxat  high  frequencies  this2recJuces2to  the  usual 
acoustic  -  gravity  dispersion  relation,  but  with  s  -*s  +  a  : 


2  2  2  2  2 
0)  (s  +  a  )  k  +  iygo)  k 
z  z 


4  2.2^  2., 

0  -  o)  (s  +  a  )k 


2  a-  M2  2o  2 

j_  +  N  s  k| 


The  vertical  component  of  the  wave  vector  is 
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There  are  acoustic  or  magnetic  pressure  waves  in  the  high  frequency  limit 
and  a  modified  gravity  wave  in  the  low  frequency  limit.  The  critical 
frequencies  are 


(A43) 


so  both  are  reduced  by  the  magnetic  field. 


For  larger  k  ,  and  low  frequencies  the  dispersion  reduces  to 


22  2  2  1  2  2  ?  9  9? 

-  uj  (s“  +  a  )  k  +  k“  s  a  k,“  +  s'k  =  0. 


This  is  the  same  as  the  magneto-gravity  dispersion  relation 

(for  horizontal  B  ) 

-  o 
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u> 


9 

k" 


2  2  2  99 

+  k  a  k_^  +  N  kj_  =  0 , 


but  with 


+  a 


Compare  fig.  14  with  Schwartz  and  Stein  (1975)  fig  3b. 
The  vertical  component  of  the  wave-vector  is 


Hence  k 

z 

critical 


2  2 
s  +a 


2  2 

2  ,  2  s  a 

^  ~  ki.  2  2 


s  +a 


■+  00  for  kx  =  a) 
frequency. 


2  2 
s  +a 

2  2 
s  a 


(A44) 


and  is  real  for  smaller  k.  and  u)  below  the 
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V.  "MODAL"  CODE 


We  wish  to  study  the  generation,  propagation,  and  dissipation  of  magneto- 
Acoustic-gravity  waves  in  order  to  determine  their  role  in  chromospheric  and 
coronal  heating,  as  well  as  to  determine  any  possible  signatures  of  different 
waves  on  line  profiles.  These  waves  are  inherently  three-dimensional  and 
three-dimensional  nonlinear  MHD  calculations  require  large  amounts  of  computer 
time.  We  are  therefore  developing  a  method  of  reducing  the  three-dimensional 
calculations  to  one  dimension.  The  essence  of  this  "modal"  method  is  to 
assume  that  the  mean  state  of  the  atmosphere  is  horizontally  uniform  then  separate 
all  the  fluid  variables  into  mean  and  fluctuating  parts,  and  expand  the  fluctuating 
parts  in  a  complete  set  of  horizontal  planforms.  The  crucial  approximation 
is  to  truncate  this  expansion  at  only  a  single  hexagonal  horizontal  planform, 
which  then  leaves  us  with  a  system  on  one-dimensional  Eulerian  partial-differential 
equations  for  the  mean  variables  and  the  amplitude  of  the  fluctuating  variables. 

For  example,  scaler  variables  are  expanded  like  the  density 
p(  K,t)  =  ~  (z,t)  +  R  (z,t)  f  (x,y) 

where  P  is  the  mean  and  Rf  the  fluctuating  part  of  the  density  and  f  is  the 
hexagonal  planform 

6 

f(x,y)  -  -  5^  e  ~i 

rT  L~ 

i=l 


(0,  |f  ), 


2JT  2  TT  -2  7T  2  IT  . 

7?S’  3S;'  -  3S  ’ 


where  S  is  the  length  of  a  side  of  the  hexagon. 
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Vector  variables  are  expanded  like  the  velocity  in  terms  of  scaler  potentials, 

fa' 1  (  fx(x,y)  U  (z,t)  +  f  (x,y)  C  (z 

V ( r , t )  =  V(Z,t)  +  a  1  <  f  (x,y)  U  (z,t)  -  f  (x,y)  C  (z,t 

y  x 

\  f(x,y)  W(z,t) 

These  hexagonal  horizontal  plan  forms  are  shown  in  Figure  19.  We  have  imposed 
this  fixed  horizontal  structure  on  the  fluid  flow.  These  1-mode  expansions 
are  substituted  in  the  equations  of  motion,  which  are  separated  into  horizontally 
averaged  and  fluctuating  parts.  For  example,  the  continuity  equation  becomes 

TT"  *  V  •  l>  *  TJ  (8W)  ■  0 

-I7  +  4-  (T  w  +  CRW)  =  a  (~  U  +  C1RU) 

3  t  3  Z 

where 

C  =  fff  =  2/ yT{~ 

C1  =  a*2  f  ff  f  +  f  fl  =3  =  1/fT 

xx  y  yJ  2  1  6 

The  advantage  of  this  method  is  that  we  have  included  the  horizontal  velocities 
and  a  crude  representation  of  the  horizontal  variation  of  the  motion,  while 
retaining  some  of  the  non-linear  -coupling  because  we  use  hexagonal  planforms, 
and  still  have  only  a  one  dimensional  system  of  equations.  The  disadvantage 
is  that  the  horizontal  strucutre  is  imposed  and  does  not  correspond  to  the 
actual  fluid  flow.  It  can,  however,  be  throught  of  as  a  standing  wave  pattern 
produced  by  six  waves  propagating  in  an  hexagonal  "pattern". 

The  modal  equations  are  solved  by  dividing  the  atmosphere  into  a  grid 
of  layers  and  transforming  the  partial  differential  equations  to  finite  difference 
equations  on  the  spatial  grid.  We  solve  the  finite  difference  equations  using 
a  two  step  Lax-Wendroff  method  (Richtmyer  and  Morton,  1969)  with  "Flux  Corrected 


Transport"  to  introduce  additional  dissipation  to  stabilize  the  scheme  (Book, 
Boris,  Hain,  1975).  We  use  the  full  one-dimensional  non-linear  Eulerian  equations 
to  test  the  solution  scheme.  At  the  moment  there  is  a  problem  with  the  FCT 
scheme  at  the  boundary  in  a  stratified  atmosphere.  We  have  tested  the  modal 
scheme  for  (a)  one-dimensional  piston  driven  shocks.  These  show  the  correct 
steepening  and  shock  formation,  but  a  bad  post  shock  profile  with  a  kink  in 
the  saw  tooth  profile.  (b)  We  are  now  calculating  waves  in  an  isothermal 
gas,  so  only  the  continuity  and  momentum  equations  need  to  be  solved,  in  order 
to  try  and  isolate  the  problem. 
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FIGURE  CAPTIONS 

Fig.  1  Chromospheric  Resonance:  Velocity  Amplitude  at  the  temperature 

_3 

minimum  for  a  given  driving  piston  amplitude,  =  3  x  10  x  sound 
speed,  as  a  function  of  driving  piston  period. 

Fig.  2  Modes:  Kinetic  energy  density  in  the  fundamental  and  harmonic 
modes  as  a  function  of  height. 

Fig.  3  Transition  region  flux  as  a  function  of  driving  piston  period 

_3 

for  a  given  driving  piston  amplitude,  =  3  x  10  x  sound  speed. 
Fig.  4  Transition  region  flux  as  a  function  of  driving  piston  amplitude 

for  driving  piston  periods  of  150  and  300  seconds. 

Fig.  5  Wave  vector  surface  at  high  frequency,  f  =  j/N  =  1.1,  in  a  strong 

magnetic  field,  a/s  =  3.16.  The  axes  are  k  ,  k  ,  k  and  full 

x  y  z 

scale  along  each  axis  is  +0.75/H,  where  H  is  the  pressure  and 
density  scale  height.  The  magnetic  field  lies  in  the  plane  of  the 
paper  at  an  angle  0  =  30l  to  the  vertical  in  this  and  all 
subsequent  3-dimensional  pictures  of  the  wave  vector  surface, 
except  figure  14.  The  spheroid  in  the  center  is  the  magnetic 
pressure  fast  mode  an  the  planes  are  the  acoustic  slow  mode 
propagating  along  the  magnetic  flux  tubes. 

Fig.  6  Wave  vector  surface  at  high  frequency,  f  =  w/N^  =  1.1,  in  a  weak 

magnetic  field,  a/s  =  0.707.  Full  scale  along  each  axis  is 
+  1.5/H.  The  spheroid  in  the  center  is  the  acoustic  fast  mode  and 
the  planes  are  the  magnetic  tension  slow  mode. 

F-g-  7  Slice  through  the  wave  vector  surface  at  high  frequency,  f  =  u)/N 

=  l.O,  in  a  strong  magnetic  field.  The  axis  are  the  real  part 

of  the  vertical,  Ke  k  ,  and  the  horizontal ,  k  .  components  of 

z  H 

the  wave  vector  -in  units  of  the  reciprocal  of  the  scale  height, 

H  The  magnetic  field  direction  is  0  =  45°  to  the  vertical 
and  f>  =  45°  to  the  plane  of  the  paper. 
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Fig.  8  (a)  Slice  through  the  wave  vector  suriace  at  high  lu-qtuikt, 

f  =  u)/N  =  1.0,  in  a  weak  magnetic  field. 

ac 

(b)  Imaginary  part  of  the  vertical  component  of  the  wave  vector, 

Lm  k  ,  along  the  same  slice  as  I  ig.  4a.  Units  arc  h  '  . 

/. 

Fig.  9  (a)  Wave  vector  suriace  ot  the  gravity  last  mode  at  low  t  requeue . , 

f  =  jj/N  =  0.5,  in  a  weak  magnetic  field, 
ac 

(b)  Wave  vector  surface  o!  the  magnetic  tension  slow  mode.  In  both 
cases  full  scale  along  each  axis  is  +  2/H. 

Fig.  10  (a)  Slice  through  wave  vector  surface  at  low  frequency,  f  =  v/N 

=  0.5,  in  a  weak  magnetic  field. 

(b)  Imaginary  part  of  along  the  same  slice.  Units  are  H 

9  and  '1  give  the  direction  of  B  . 

o 

Fig.  11  (a)  Slice  through  wave  vector  surface  at  low  frequency,  f  =  .../N 

=0.5,  in  a  nearly  vertical,  9  =  3",  magnetic  field. 

(b)  imaginary  part  of  k  along  the  same  slice.  Units  are  H 

Fig.  12  (a)  Slice  through  wave  vector  surface  at  low  frequency,  f  =  w/N 

c 

=  0.5,  in  a  weak  magnetic  field  perpendicular  to  the  B  -  g  plane, 

o 

■t>  =  90  . 

(b)  imaginary  part  of  k  along  the  same  slice.  Units  are  H  . 

7, 

Fig.  13  (a)  Wave  vector  surface  of  "slow  gravity  mode"  at  low  frequency, 

f  =  /N'  =  0.5,  in  a  moderate  magnetic  field, 

ac 

(b)  Wave  vector  surface  of  magnetic  "fast"  mode  in  same  case.  Full 


scale  along  each  axis  is  +  1/U. 
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Fig  .  Ill 


Fig.  1b 


Fie.  16 


Fig.  17 


Fig.  18 


7ig.  19 


Wave  vector  surface  at  low  f  ro<|uencv ,  I  =  . /N  =  0.5,  in  a  strong 
magnetic  field.  Bulge  is  magnetic  pressure  last  mode.  Other  mode 
is  a  plane  through  origin.  Full  scale  along  each  axis  is  +  O.l/H. 

(a)  Slice  through  wave  vector  surface  at  low  frequency, 

f  =  o)/N  =0.5,  in  a  strong  field, 
ac 

(b)  Imaginary  part  of  k?  along  the  same  slice.  Units  are  H  h 

(a)  Slice  through  wave  vector  surface  at  low  frequency, 

f  =  rn/N  =  0.5,  in  a  strong  magnetic  field,  in  a  direction 
ac 

perpendicular  to  the  B^-g  plane,  4>  =  90°. 

(b)  Imaginary  part  of  along  the  same  slide.  Units  are  H 

(a)  Slice  through  the  wave  vector  surface  at  low  frequency 
f  =  „  =  0,5»  in  3  nearly  vertical,  9  =  5°,  strong  magnetic 

clC 

field. 

(b)  Imaginary  part  of  k  along  the  same  slice.  Units  are  H_1. 

z 

Wave  vector  surface  for  nearly  horizontal  magnetic  field,  0=  85°, 

at  frequency,  f  =  w/N  =  0.5,  in  a  moderate  magnetic  field. 

ac 

2  2 

cutoff  frequency  in  this  case  is  N  s/(s  +  a  ) 2  =  0.^8  N 

ac  ac. 

The  inner  fast  mode  surface  is  similar  to  the  magneto-gravity 
mode  surface  •  (Schwartz  and  Stein,  1975  fig-  3b). 

Hexagonal  planform  f(x,y). 
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FLUX  vs.  PISTON  AMPLITUDE 


Figure  4 


Figure  5 
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Figure  17b 
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